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Abstract 

We study the effect of chemotactic signahng among mesenchymal cells. We show that the 
particular physiology of the mesenchymal cells allows one-dimensional collapse in contrast to the 
case of bacteria, and that the mesenchymal morphogenesis represents thus a more complex type 
of pattern formation than those found in bacterial colonies. We finally compare our theoretical 
predictions with recent in vitro experiments. 
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The development of spatial patterns is one of the most important topics in embryology. 
The formation of structure in embryology is known as morphogenesis. Although genes play 
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a crucial role in the control of pattern formation, the importance of the mechanochemica 
teractions among the cells and their environtment has been recognised in several works 
One of the advantages of this approach is that it has the potential for self correction in con- 
trast to the Turing chemical prepattern approach. Embryonic development is usually a very 
stable process with the embryo capable of adjusting to many outside disturbances. The 
prepattern approach implies the existence of potentially unstable processes and make it dif- 
ficult for the embryo to make the necessary adjustement to such disturbances as development 
proceeds 

In this work we are concerned with one type of early embrionyc cells known as dermal 
or mesenchymal cells, responsible of the formation of highly organised patterns on skin 
such as the primordia which become feathers and scales, and the condensation of cells 
which mirror the cartilage pattern in developing limbs. Mesenchymal cells are capable of 
independent movement, due to long finger- like protusions called filodopia which grab onto 
adhesive sites and pull themselves along: spatial aggregation patterns in these appear as 
spatial variations in cell number density 0,0]. These cells can also secrete fibrous material 
which helps to make up the extracellular matrix tissue within which the cells move. However, 
experimental evidence indicates that there is not such secretion during chondrogenesis and 
pattern formation of skin organ primordia j^, so we will neglect this contribution to the 
dynamics. 

Here, we will analize the role of chemotaxis in mesenchymal morphogenesis. It is known 
that chemotactic signaling is one of the most important mechanisms that lead to pattern 
formation in bacterial colonies Q], suggesting that its role might be crucial in morphogene- 
sis. Actually, the presence of a powerful chemoattractant has been identified as one of the 
active responsibles of pattern formation in mesenchymal self-organization 3|- Probably, the 
simplest mathematical model for chemotactic aggregation is the Keller-Segel model 

dtp = DhV^p-V-ikpVc), (1) 
dtc = D^V^c + ap. (2) 

Here Df, is the cellular diffusion constant, k the chemotactic coefficient, a the rate of attrac- 
tant production, and Dc the chemical diffusion constant. The terms in Eq.(^ include the 



diffusion of the cells and chemotactic drift. Eq.Q expresses the diffusion and production of 
attractant. Nondimensionalizing system ()1I2|) we get: 

dtp = V^p-V -{pVc), (3) 
edtC = V^c + p, (4) 

where e = Dt/Dc- An efficient chemotactic communication implies that the diffusion of the 
cells is much slower than attractant diffusion, which leads to consider e = 0. We finally 
arrive at the following nonlinear partial differential system: 

dtp = V'p-V-ipVc), (5) 
- V'c = p - fco, (6) 

where = ^ pdx, and Q is the region of the space where the system is defined, \Q\ being 
its volume. Note that the introduction of ko provides a solvability condition for Eqs. ()5|6|l 
in the case of no ffux boundary conditions jsf. This system is known to blow up in finite 
time for dimension d> 2, but all the solutions are regular for = 1 This means that 
in a three-dimensional system, while collapse to infinite density lines and points can occur, 
collapse to an infinite densi ty s heet is mathematically impossible. This fact crucially affects 



the patterns that can form Il2j|. Actually, both types of chemotactic collapse have been 
already observed in experiments performed with Escherichia Coli 3, • 

In the case of mesenchymal cells, far more complex than a bacteria like Escherichia 
Coli, the situation gets more involved. The supposition of short range diffusion (or simply 
diffusion), that applies well to dilute systems, it is not, in general, sufficiently accurate in 
such systems in which the cell densities are relatively high. The long filopodia extended 
by the cells can sense density variations beyond their nearest neighbours and so we must 
include a nonlocal effect on diffusive dispersal since the cells sense more distant densities 
and so respond to neighbouring averages as well [3]. 

Long-range diffusion was tradionally modeled by the inclusion of a biharmonic term of 

the form V^. This comes from the known fact that: 

2 (p(x,t)) -p(x,t) 
V p oc — , as R^O, (7) 

where (p) is the average density in a sphere of radius R about x, that is 

(p(x,t)) = [ p(x + r,t)c^r, (8) 



where V is the sphere of radius R. Because R ^ 0, this suggests in the one-dimensional 
case the following Taylor series expansion: 



p{x). 



r ( l + -rX + - I d: 



(9) 



In the case of diffusion in an isotropic medium, after integration and truncation after the 
fourth term we obtain: 



(p(x,t)) -p(x,t) 

i?2 



{D2dl + R^Did^p + o{R^ 



(10) 



where the average is performed over the closed interval [— -R, R] . This way we get the 
following extended diffusion equation: 



dtp = {D^dl + R^D^dDp, 



(11) 



with D2,D4 > 0. An initial value problem to Eq.()ll|) blows up at finite time; this is a 
consequence of the assymptotic character of the functional Taylor expansion ()10|) . Physically, 
this means that the cells move randomly but up a cell density gradient, a fact that goes 
against experiment. Further, the next higher truncation leads to a better behaved equation, 
but the corresponding solution is negative somewhere as known from Pawula's work jlfil ]. 
Also, any aproximation beyond the second order leads to a nonphysical increase of the 
number of boundary conditions, so we have to conclude that this is not a proper way to 
generalize diffusion. We can solve this problem regrouping the terms in Eq. lfTIH) in the 



manner of Fade to get 
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id: 



{p{x,t)) - p{x,t) 



R^ l-R'^{Di/D2 
The resulting extended diffusion equation is then: 



i? 0. 



(12) 



dtp 



l-R\D^/D2) 

where the diffusive operator is to be interpreted in the Fourier transform sense 



(13) 



Eq. (fT^ seems to be a proper extension of the diffusive aproximation as shown in Ref . [3| • 

We will see that considering long-range diffusion in system ()5I6|) will lead to a finite time 
singularity in d = 1, that implies collapse to an infinite density sheet in a three-dimensional 
system and to an infinite density line in a two-dimensional system. We are going to show this 
fact analitically, since numerical calculations are extremely unstable to precisely compute 



the existence of blow-ups in partial differential equations 
We will thus consider the system: 

V2 



= 1 -'e2V2 ^ - Vc ■ Vp + - kop, 

- V'c = p-ko, (15) 

in one spatial dimension. Here, e is proportional to the mean radius of a cell and the natural 
boundary conditions are no flux boundary conditions (with a long-range gradient in the case 
ofp) 

0, (16) 



X 



an 

where Q is the closed interval Q = [—L,L], dQ being its boundary. Note that integrating 
the equation 

dtp = YZ^2P - ■ (P^-c) (17) 
over VL and applying the boundary conditions we get the conservation of the total mass 
of p (as it should be since we are only considering movement of the cells) and thus the 
conservation in time of k^. To clarify the notation, let us explicity write the norm of a 
function / belonging to a L^'iVL) space, 1 < p < oo: 

i/p 



Lp{n) 



IfYdx] . (18) 



From equation ()T5j) we get: 



— / pdxcdxpdx + / p^dx — ko p^dx. (19) 
Jo. Jn Jq 

Now, we are going to estimate all the terms appearing in the right hand side of this equation. 
Integrating by parts the second term in the right hand side of Eq. lfTIH) : 

pd^cdxpdx = P^f^xc|g^ - / dxpdxCpdx 



- I pdlcpdx, (20) 
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that implies 



pdxcdxpdx = ^7: P d^cdx 
n ^ Jq 



p^dx. 



The first term in the right hand side of Eq. ()19|) can be estimated as follows: 



dl 



pdx < 



dl 



2Q2 



pdx 



< 



P- 



1 - e^dl 



dx < \ \p\\^2 



^x 



P 



L2(n) 



(21) 



(22) 



where we have used Holder's inequality (see below). By performing the shift of variables 
y = x/e, we get: 



^ ^x 



.(3/2) 



y 



P 



< 



L2(n/e) 



- ^2 

y 
N 

\ \P\\L^{Q/e) ' 



(23) 



dlil — dl) ^|. Let us clarify a bit this last step. We have used the fact that the 



y^ y> 

2(^ _ y2-j-i -g bounded on every space, with I < p < 00. This means that 



where 

operator V^(l 

we can assure that ||V^(1 — ^'^)^^ f\\LP{n) — ^ \ \f\\Lp(n) ^'^^ every / belonging to Lp{Q) and 
a constant that does not depend on / (and thus is called the norm of the operator). 
This fact can be easily seen once one realizes that the Fourier transform of the operator 
V2(l — V^)"^ is a bounded function of the wavevector, and a rigorous proof can be found 

n 

m m. We can again shift variables x = ey to get: 



N 



P 



1 - 5292 



pdx < ||p||i2(o) • 



Finally, we can conclude our estimate as follows: 



92 



^2Q2 

X 



pdx > 



92 



2Q2 



pdx 



> 



N 
'I2 



|2 

Il2(Q) 



(24) 



(25) 



Now we are going to estimate the third term in Eq. p9j) : 



37 _ 1 1 1 13 

p ax — \ \P\\l3(^q) ■ 



(26) 
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Holder's inequality reads (for a rigorous proof of Holder's inequality see j2d|) 



\uv\ dx < lbllL9(n) ' 



l<P,g<oo, - + - = 1. 27) 
P g 

Choosing t> = 1 we get: 

/ |m| < C , (28) 

where C = \VL\^^'^ . With this estimate we can claim that: 



^ - ' p^dx <C\\p^ 



(i/p) 



c(^j p^^dx 



/ r \(2/3) 

Ci^j^p'dxj =C||p||j3(^), (29) 
where we have chosen p = 3/2 (and correspondingly q = 3). This implies that: 

\\p\\LHn)> DMIl^q)^ (30) 
where D = Therefore, we have the final estimate: 

^||p|li2(f^) > ^ (||p|lL2(f^))^ ^ ^ - B \\p\\l,^^^ , (31) 

where A, B > are constants, A = l^l]^^/"^ and B = ^ + k^. We are thus going to study 
the dynamical system: 

^ = Ax'/^ - Bx. (32) 

dt ^ ' 

This system has two fixed points, x = and x = (B/A)^ > 0. A linear stability analysis 
reveals that the positive fixed point is linearly unstable, meaning that every initial condition 
Xq > {B/AY will stay above this value for all times. Further, we know that the solution will 
grow without bound in this case, so we can claim the existence of two constants, Iq < oo 
and < Co < A, such that Ax'^/'^if) - Bx{t) > C^x^/'^it) for every t > to. This implies that 



^ Ilp|li2(f^) > Co (||p|li2(f^))^ ^ ^ (33) 



for t > to, and for an adecuate initial condition. Solving this equation gives: 

\\pi-Ml^n) > , \ — (34) 



for t > ti > to, and for an adecuate initial condition. And every adecuate initial condition 
must fuUfill 

\\pi-M\lHn)>^m 
+ ^IIP(-,0)||^.(^)+ l^||p(-,0)||i.(^), (35) 

like, for instance, p{x,0) = (x^ + 5)"^/^ and S small enough. Thus we are finally led to 
conclude that the system does blow up in finite time. 

It is interesting to see that a higher order correction to the long-range diffusion does not 
alter this behaviour. Indeed, the next order corresponds to an operator of the form 

^2 



(36) 



i-cai 

The Fourier transform of this operator is bounded in the long wavelength limit, meaning 
that we will get an analogous result for this case, indeed 

(_l.2\ 1 + ^ 

l-C(-fc^) 

Actually, further generalizations of it, containing an odd number of fractions between an 

even number of expressions that have a second order derivative will have a well-defined long 

wavelength limit, and will thus lead to the same result. 

It has been argued that each of the aggregates in a p attern corresponds to a density 

singularity in the hydrodynamic description of the cells [uj . Our analysis predicts a different 

way of pattern formation from the usual models of chemotactic aggregation. In particular, 

an initial diffusive band can form a singularity by collapsing only one of its dimensions to 

zero thickness. This type of one-dimensional collapse has already been empirically observed: 

in vitro experiments showed that mesenchymal cells are able to aggregate by collapsing only 

n 

one of the dimensions of the culture into stripes [7] . After a few days, an initial homogeneous 
layer begins to develop spatial structure, the cells beginning to align with their neighbors 
to form "swirls". This might be related to the fact that the homogeneous layer is not the 
" adecuate" initial condition that we derived in our theoretical analysis. They can in contrast 
aggregate by collapsing two spatial dimensions following a standard Keller-Segel mechanism. 
When the distribution of cells is driven far enough from the homogeneous distribution, the 
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"adecuate" initial condition is then achieved, and the culture of cells begins to aggregate by 
collapsing only one of its dimensions into ridges. 
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